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We demonstrate a degenerate parametric oscillator in a silicon-nitride microresonator. We use two frequency- 
detuned pump waves to perform parametric four-wave mixing and operate in the normal group-velocity- 
dispersion regime to produce signal and idler fields that are frequency degenerate. Our theoretical modeling 
shows that this regime enables generation of bimodal phase states, analogous to the y^^^-based degenerate 
OPO. Our system offers potential for realization of CMOS-chip-based coherent optical computing and an 
all-optical quantum random number generator. 
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Researchers recently proposed using degenerate op¬ 
tical parametric oscillators (OPO’s) to realize a novel 
form of coherent optical computing by simulating the 
quantum Ising model m- Such an approach offers 
the possibility of solving non-deterministic polynomial 
time (NP)-hard computation problems [1]. A first proof- 
of-principle demonstration based on a degenerate 2nd- 
order nonlinearity OPO using bulk optics offers the 
promise that this new paradigm could be realized ISIS]. 
The approach is based on the non-equilibrium phase 
transition that occurs at threshold for optical parametric 
oscillation, in which the generated held with respect to 
the pump held locks to one of two possible states offset 
by TT. Below threshold, the signal held is described by 
a squeezed vacuum state, where the strength of squeez¬ 
ing depends on the ratio of the intracavity power to the 
threshold power 00. The vacuum huctuations of the 
out-of-phase quadrature are suppressed, while those of 
the in-phase quadrature are enhanced. It is along the 
in-phase quadrature that, at threshold, the system bi¬ 
furcates into two possible states with opposite phases. 
By utilizing a network of coupled OPO’s, more com¬ 
plex phase-locked states can be achieved which encode 
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the ground-state of the quantum Ising model [T] , where 
the coupling strengths of its spins are equivalent to the 
coupling between the OPO’s. It has been suggested 
that such an OPO system offers potential for a sim¬ 
ulation device that can outperform any other classical 
computation [5]. Recently, a 4-OPO system based on 
the nonlinearity has been demonstrated [2]. Alter¬ 
natively, the bi-phase state output can be used to realize 
a quantum random number generator that compared 
to other implementations requires only classical signal 
detection (no single photon detection), is intrinsically 
unbiased, and thus does not need any post-procession 
of the raw output, which has applications in cryptogra¬ 
phy 0, Monte Carlo simulations and quantitative 
finance m- The scalability and robustness of these cou¬ 
pled OPO’s will be crucial for these applications. 

An alternative approach suggested in [5] is to utilize 
a nonlinearity to realize an analogous degenerate 
OPO with dynamics that would mirror the sys¬ 
tem. The silicon-based photonics platform appears to 
be highly suited for this application since it allows for 
a robust, compact, integrated system that offers feasi¬ 
bility for scaling and stability. In recent years, there 
have been various nonlinear optics demonstrations in 
microresonators [HHizi. In particular, there has been 
signihcant development in frequency comb generation 
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Fig. 1. (a) Four-wave mixing energy level diagram in which 
two pump pumps are converted to two frequency-degenerate 
photons, (b) Scheme for our OPO in Si 3 N 4 microres¬ 
onator, where Ai and A 2 are the pump fields, As and Ai 
are the generated signal and idler. 5 is the frequency offset 
between the pumps and the center frequency, A is the fre¬ 
quency offset from degeneracy between the signal/idler and 
the center frequency. 
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through optical parametric oscillation in high-Q mi¬ 
croresonators [m [Mni [2nH26] that utilizes paramet¬ 
ric four-wave mixing (FWM) interactions based on the 
third-order nonlinearity The silicon nitride (Si 3 N 4 ) 
platform offers potential as a platform for coupled de¬ 
generate OPO’s due to its CMOS compatibility, low 
linear and nonlinear losses at telecommunication wave¬ 
lengths, and high nonlinearity and dispersion engineer¬ 
ing [25]. Recently, there have been demonstrations of 
inverse FWM using optical fibers in the normal group- 
velocity dispersion (GVD) regime [23], and Inagaki, et 
al. [30] has shown in initial experiments that a dual- 
pumped fiber OPO’s multiplexed in the time do¬ 
main can be used to simulated a 1-D Ising model. 

In this paper, we theoretically and experimentally 
demonstrate that a dual-pumped, microresonator-based 
x(3) OPO can achieve degenerate oscillation. In our 
scheme, a Si 3 N 4 microresonator is dual-pumped to 
generate a frequency-degenerate signal and idler pair 
through parametric FWM [Fig. [^a)]. We perform theo¬ 
retical analysis to determine the conditions under which 
degenerate oscillation and the associated random binary 
phase state is achieved and find that it requires normal 
GVD using a dual-pump. We show experimentally that 
frequency-degenerate oscillation is achieved by control¬ 
ling the frequency offset between the two pumps. 

Our configuration consists of two pump fields denoted 
by Al and A 2 that are offset by frequency ±<5 from the 
degeneracy point [see Fig. ^b)]. There has been a 
number of experimental studies that have investigated 
parametric FWM with two pump waves [STIIS^ . En¬ 
ergy conservation implies that the generated signal and 
idler fields with amplitudes denoted by Ag and re¬ 
spectively, can grow at a frequency offset A from degen¬ 
eracy for the signal/idler pair. In order for degenerate 
oscillation to occur, the parametric gain for the signal 
and idler must peak at A = 0. We model the nonlinear 
interaction using the nonlinear Schrodinger equation. 
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Fig. 2. Parametric gain profile assuming two undepleted non¬ 
degenerate pumps. Dashed lines indicate the position of the 
pumps. Pump separation is 14 nm and the center wavelength 
between the pumps is 1550 nm. (a)-(c) show the gain profiles 
in normal GVD {P 2 = 193 ps^/km) for the ratio between Lo 
and Lml of 2, 0.5 and 0.2, respectively. The nonlinear co¬ 
efficient (7 = 0.94 W“^m“'^) is based on the 690x1300 nm 
cross-section Si 3 N 4 waveguide used in our experiment, (d)- 
(f) show the gain profiles in anomalous GVD {P 2 = —193 
ps^/km) for the same ratios. The ratio between Ld and 
Lnl is adjusted by changing the pump power. The effect 
of higher-order dispersion terms is negligible because of the 
small frequency separation between the pumps as well as the 
waveguide dispersion. 


where /?2 is the GVD parameter and 7 is the nonlinear 
parameter proportional to x^^^- We write the total field 
with respect to the center mode frequency as 


A(z,t) = Ai{z)e^^'^+A2{z)e-^^^+Asiz)e^^^+A,iz)e-^^\ 

( 2 ) 

We assume that the pump waves are undepleted and 
have equal powers P = |Ai|^ = |A 2 |^. The solution 
for the signal/idler pair has an exponential form 
where g is the gain parameter. 


g = i/lG-fP - ^ 2 {S^ - A^)] X [ 27 P + ^ 2 (< 52 -A2)]. 

(3) 

Our analysis indicates that with suitable control of the 
pump power and offset frequency, it is possible to find a 
condition for achieving parametric gain at modest nor¬ 
mal GVD (i.e., f32 >0), which can be readily achieved by 
engineering the waveguide cross-section m- We inves¬ 
tigate the parametric gain profile based on the param¬ 
eters of a 690x1300 nm cross section Si 3 N 4 waveguide 
pumped with two single-frequency lasers. Figure[^shows 
the dependence of the gain on the dispersion length 
Ld = 1 / 5 ^ 1 / 32 ! and the nonlinear length LpfL = 

We show three different regimes, each with a different 
ratio between Ld and Ljv l for both normal and anoma¬ 
lous GVD and find that maximum parametric gain can 
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be achieved at the degenerate signal/idler frequency only 
when Ld > Lml in the normal GVD regime [Fig. [^a)]. 
Operating in this well-defined normal GVD regime is 
essential in order to implement the bi-phase state, as 
this ensures that the degenerate pairs have the lowest 
threshold for oscillation. 

Furthermore, in order to ensure degenerate oscillation 
in a microresonator cavity, the position of the two pumps 
and the signal/idler within the cavity resonance is criti¬ 
cal. Figure shows the resonator conditions in the nor¬ 
mal GVD regime. We define the resonator mode fre¬ 
quencies near the frequencies of pump 1, pump 2, and 
the degenerate signal/idler frequencies as /i, / 2 , and 
/r, respectively (indicated in red in Fig. [^. The cen¬ 
ter frequency between the two pump resonance modes 
/c = is offset from the signal/idler resonance mode 

by A/ = fr — fc- This frequency offset is dependent on 
the dispersion and, for a cold cavity with normal GVD, 
it is given by A/coid = 'Cg/dTrL^i, where Vg is the group 
velocity. In a hot cavity, the frequency offset is due to 
both dispersion and the Kerr nonlinearity and is given 
by A/hot = - lTf)- frequencies of 

/pp and /pp (indicated in blue in Fig. [^, we define 
detuning of the pump frequency from the peak of the 
resonance as Sf = fp,i — fi = fp.2 — h- The signal/idler 
position with respect to the resonance for the zero de¬ 
tuning case, where both pumps are tuned to the peak 
of the respective resonances, is shown in [Fig. da)]. As 
discussed earlier, for the degenerate signal/idler to os¬ 
cillate, Ljj > Lml must be satished, which corresponds 
to a negative A/hot that could jeopardize the oscillation 
of the degenerate pair. If [A/hot! > Si'/2, where is 
the cavity linewidth, the degenerate pair is suppressed. 
However, we can take advantage of pump detuning to 
achieve the optimum oscillation conditions in the res¬ 
onator. First, since Ljsil is dependent on the total power 
built-up inside the resonator, pump detuning affects the 
parametric gain profile required for oscillation. In addi¬ 
tion, the pump frequency tuning can be used to control 
the optimal signal/idler frequency position provided that 
I A/hot — Sf \ < (5i//2 [Fig. I^b)]. This allows for all fields 
to lie within their respective cavity resonances allowing 
for degenerate oscillation to occur. In practice, the ulti¬ 
mate detuning range that can be achieved is limited by 
the thermal effect in the resonator from pumping two 
different resonances. 

We numerically simulate the degenerate OPO using 
the modified Lugiato-Lefever equation [35H38], which 
we previously developed for modeling parametric fre¬ 
quency comb generation in microresonators |39j . The 
model includes effects of higher-order dispersion and 
self-steepening, which enables simulations over a wide 
range of conditions. The dispersion of the microres¬ 
onator mode is simulated using a finite-element mode 
solver [Fig. |^a)]. The simulated spectra of the OPO is 
shown in Fig. ®b) for a Si 3 N 4 microresonator with a 
free spectral range (FSR) of 200 GHz and a cross sec¬ 
tion of 690 X1300 nm, which corresponds to normal GVD 
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Fig. 3. Resonator condition for dual-pump OPO. The rel¬ 
ative positions of the resonances are shown for (a) a hot 
cavity with zero pump detuning and (b) hot cavity with 
pump detuning, /i, / 2 , and fr (red) are the resonance fre¬ 
quencies for pump 1, pump 2, and degenerate signal/idler 
modes, respectively. /p,i, /p,2, fs,i (blue) are the two pump 
frequencies and the generated signal/idler frequency, respec¬ 
tively. fc = A/hot = fr — fc, and the pump detuning 

= /p,i - /i = /p,2 - /2- 


(a) 


a 

> 

CD 


(b) 


(c) 


o 

CL 


400 



-200 


1200 1300 


1400 1500 1600 

Wavelength (nm) 


1700 1800 



1450 


1500 1550 1600 

Wavelength (nm) 


1650 


■“t 


10 


20 


30 40 

Time (ns) 


50 


60 


70 


Fig. 4. (a) Simulated GVD for the fundamental TM mode 
of 690x1300 nm cross section Si 3 N 4 microresonator. Inset: 
Si 3 N 4 microresonators, (b) Simulated intracavity spectrum 
for degenerate OPO in Si 3 N 4 microresonator. Inset shows the 
spectrum on a dB scale. The pump offset from the degenerate 
frequency S = 4xFSR. (c) Evolution of the phase of the 
signal for two realizations. Bi-phase states are achieved with 
a TT phase offset. 


at the pump wavelengths of 1535 and 1565 nm for the 
fundamental TM mode. The two pumps are each offset 
from the frequency degeneracy point by 5 = dxFSR’s. 
The simulated spectra show that parametric oscillation 
of the degenerate signal/idler pair occurs only under 
pump detuning conditions consistent with our analysis 
above. Multiple simulations with different initial noise 
conditions indicate that the degenerate signal held set¬ 
tles into one of two phase states which are offset by tt 
[F ig. i c)], which is precisely the required behavior for 
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Fig. 5. OPO spectra from two pumps for pump offsets S of 
(a) 6.5XFSR, (b) 5xFSR, and (c) 4xFSR. The FSR is 200 
GHz. The pump offset of 4xFSR (c) enables the generation 
of a degenerate signal/idler pair between the pumps, which 
results in the generation of the bi-phase state. 


simulating the Ising model as well as for quantum ran¬ 
dom number generation. Since the generation rate is 
largely dependent in the cavity lifetime, the current sys¬ 
tem should be capable of producing random numbers at 
rates of 100 MHz. We believe that by reducing the size 
of the resonator and with suitable tuning of the cavity 
linewidth, random number rates approaching 1 GHz can 
be readily achieved. 

In our experiments, we use a Si 3 N 4 microresonator 
with an effective resonator cross section of 690 x 1300 nm, 
which allows for the necessary normal GVD for the fun¬ 
damental TM polarization mode. The microresonator 
has a loaded Q of 3 million and an FSR of 200 GHz. We 
combine the output beams from two amplified tunable 
single-frequency pump lasers using a 50/50 fiber coupler 
and inject them into the Si 3 N 4 bus waveguide which is 
coupled to the microresonator. WDM filters with 9-nm 
bandwidth are used to remove the amplified spontaneous 
emission from the erbium-doped fiber amplifier at each 
of the respective outputs before combining. The input 
polarization is set to quasi-TM. The output from the 
device is collected using a free-space objective and sent 
into a 80/20 fiber coupler such that 20% of the signal 
is sent to an optical spectrum analyzer for characteri¬ 
zation, and the remaining light is sent to two different 
WDM filters and a photodiode to monitor independently 
the transmission of each of the pumps to ensure that 
they are coupled to a cavity resonance. The combined 
pump power in the coupling waveguide is 125 mW. 

In order to generate the dual-pumped OPO, we tune 
each pump into a microresonator cavity resonance while 
taking into account the thermal frequency shift that oc¬ 
curs as the pump power builds up in the cavity. The 
reduction of the pump transmission along with cascaded 
FWM due to the contribution from the two pumps in¬ 
dicates that the pump is coupled to a cavity resonance. 


Figure shows the generated OPO spectra for various 
pump frequency offsets S from the degeneracy point of, 
6.5xFSR [Fig. [ga)], 5xFSR [Fig. [gb)], and 4xFSR 
[Fig. i c)], which corresponds to pump separations of 
13x, lOx, and 8xFSR, respectively. Each of the gener¬ 
ated spectra shows a gain peak between the two pump 
wavelengths, consistent with the theoretical results. As 
discussed earlier, the pump separation along with the 
pump detuning within the resonance must be carefully 
chosen to achieve conditions for degenerate oscillation. 
Figures a) and (b) show multiple spectral lines gener¬ 
ated between the two pumps, with spacings of lx FSR 
and 2xFSR, respectively. This spacing depends on the 
first lines that oscillate near the frequency degeneracy 
point. Thus, if the degenerate signal/idler is the first 
spectral line that reaches the oscillation threshold, there 
will be no other lines generated between the pumps. For 
the offset of <5 = 4xFSR, which corresponds to pump 
wavelengths of 1545 nm and 1557.7 nm [Fig. i c)], we 
generate a single frequency component between the two 
pumps, corresponding to a degenerate signal/idler pair 
at the center frequency, which will enable the generation 
of the bi-phase state. 

In conclusion, we theoretically and experimentally 
demonstrate frequency-degenerate parametric oscilla¬ 
tion via dual-pump FWM in Si 3 N 4 microresonators 
operating in the normal GVD regime. We show that 
the Si 3 N 4 platform offers the potential of developing 
a scalable GMOS-compatible platform that utilizes 
the unique properties of light to solve challenging 
computational problems and for use as a system for 
random number generation. 
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